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ABSTRACT: The main objective of the work is to analyse the performance and flow characteristics of convergent 
divergent nozzle and also to compare the numerical values of the two methods i.e "HIT & TRIAL METHOD" AND 
"ANALYTICAL METHOD". In this paper we have determine the location and strength of normal shock wave in the 
divergent portion of the nozzle under varying operating conditions and with different nozzle geometry. 
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I. INTRODUCTION 

A nozzle is a relatively simple device, just a specially shaped tube through which hot gases flow. However, the 
mathematics, which describes the operation of the nozzle, takes some careful thought. Nozzles come in a variety of shapes 
and sizes. Simple turbojets, and turboprops, often have a fixed geometry convergent nozzle as shown on the left of the 
figure. Turbofan engines often employ a co-annular nozzle. The core flow exits the centre nozzle while the fan flow exits the 
annular nozzle. Mixing of the two flows provides some thrust enhancement and these nozzles also tend to be quieter than 
convergent nozzles. Afterburning turbojets and turbofans require a variable geometry convergent -divergent - CD nozzle. In 
this nozzle, the flow first converges down to the minimum area or throat, then is expanded through the divergent section to 
the exit at the right. The variable geometry causes these nozzles to be heavier than a fixed geometry nozzle, but variable 
geometry provides efficient engine operation over a wider airflow range than a simple fixed nozzle. Rocket engines also use 
nozzles to accelerate hot exhaust to produce thrust. Rocket engines usually have a fixed geometry CD nozzle with a much 
larger divergent section than is required for a gas turbine. 

1.1 SHOCKS IN NOZZLE 

A shock wave (also called shock front or simply "shock") is a type of propagating disturbance. Like an ordinary 
wave, it carries energy and can propagate through a medium (solid, liquid or gas) or in some cases in the absence of a 
material medium, through a field such as the electromagnetic field. Shock waves are characterized by an abrupt, nearly 
discontinuous change in the characteristics of the medium. Across a shock there is always an extremely rapid rise in 
pressure, temperature and density of the flow. In supersonic flows, expansion is achieved through an expansion fan. A shock 
wave travels through most media at a higher speed than an ordinary wave. Unlike solutions (another kind of nonlinear 
wave), the energy of a shock wave dissipates relatively quickly with distance. Also, the accompanying expansion wave 
approaches and eventually merges with the shock wave, partially cancelling it out. Thus the sonic boom associated with the 
passage of a supersonic aircraft is the sound wave resulting from the degradation and merging of the shock wave and the 
expansion wave produced by the aircraft. When a shock wave passes through matter, the total energy is preserved but the 
energy which can be extracted as work decreases and entropy increases. This, for example, creates additional drag force on 
aircraft with shocks. Shock waves can be 

□ Normal: at 90° (perpendicular) to the shock medium's flow directions. 

□ Oblique: at an angle to the direction of flow. 

□ Bow: Occurs upstream of the front (bow) of a blunt object when the upstream velocity exceeds Mach 1. 

Shock waves form when the speed of a gas changes by more than the speed of sound. At the region where this occurs sound 
waves travelling against the flow reach a point where they cannot travel any further upstream and the pressure progressively 
builds in that region, and a high pressure shock wave rapidly forms. Shock waves are not conventional sound waves; a shock 
wave takes the form of a very sharp change in the gas properties on the order of a few mean free paths (roughly micro- 
meters at atmospheric conditions) in thickness. Shock waves in air are heard as a loud "crack" or "snap" noise. Over longer 
distances a shock wave can change from a nonlinear wave into a linear wave, degenerating into a conventional sound wave 
as it heats the air and loses energy. The sound wave is heard as the familiar "thud" or "thump" of a sonic boom, commonly 
created by the supersonic flight of aircraft. 

II. METHODOLOGY 

2.1 ANALYTICAL APPROACH 

To determine shock location and shock strength in convergent divergent nozzle. 

STEP 1: Find the pressure ratio that will produce a shock in the divergent portion of the nozzle. 
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STEP 2: Determine exit Mach number (M e ) 
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STEP 3: Use M e to determine P e /P 0 
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STEP 4: Since M e <l, P e =P b 
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STEP 5: Determine M b by using the value of P02/P0 
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STEP 6: Determine Shock Location A s /A t 
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STEP 7: Determine Shock Strength 
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STEP 8: Determine Temperature Ratio across the shock 
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STEP 9: Determine Pressure Ratio across the shock 
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STEP 10: Determine density ratio across the shock 
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2.2 HIT AND TRIAL APPROACH 

To determine the shock location and shock strength in convergent divergent nozzle 

Let us consider a convergent divergent nozzle with inlet and outlet section specified in the diagram as 1 and 5 
respectively. In the diagram section 2 represent the throat while the section 3 and 4 represent the flow condition before and 
after the shock. 
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Fig. 2. 2.1 Convergent divergent nozzle 



Let, 

A e /A t =exit area/throat area=1.53 
Inlet condition P 0 =latm. 

Now let us determine the various flow characteristics for this case. 



STEP 1: 

If we assume chocked flow at the throat then from the IFT for A e /A t we have two solutions 
1) For FIRST critical point(fully isentropic subsonic flow) 

From IFT at (A e /A t )=1.53,we have 

(P e /P 0 )=.886 

M e =.42 

For THIRD critical point(fully isentropic supersonic flow) 
From IFT at (A e /A t )=1.53,we have 
(P e /P 0 )=.154 
M e =1.88 

Now as we know that normal shock takes place in the divergent portion of the nozzle i.e in the supersonic flow. Therefore 

for supersonic flow conditions we have 

M e =1.88 and 

P e /P 0 (P 3 /Po3)=.154 

STEP 2: 

Now from the normal shock table for M e =1.88, we have 

M e (aftershock)=.599 

P y /P x (after shock)=3.957 

And the operating pressure ratio will be 

Prec/P 0 =P4^01=P4^3*P3/P03*P03/P01 

P 4 /P 01 =3.957*.154*l 
P4/P 0 i=.609378 

Thus for our C-D nozzle with A e /A t =1.53,any operating pressure ratio between .886 and .609378 will cause a normal shock 
to be located some where in the divergent portion of the nozzle. 
STEP 3: 

Now let us find shock location at the operating pressure ratio of 0.75 
Let us assume shock wave is located at A s /A t = 1.204 

Note-this value is selected because it is one of the numbers in IFT. by selecting this value we don't need to interpolate. 

Key equation: 

P e =P e /P 02 *P 02 /P 01 *P 01 

From IFT, corresponding to A s /A t = 1.204, we have 
Mi=1.54 

From NST corresponding to M^l.54, we have 

M 2 =.687 and P 02 /Poi=.917 

From IFT corresponding to M 2 =.687,we have 

A s /A*=1.1018 

.A.g/.A. — .A.g/.A.£ * .A.^/ .A. § * .A. s/.A. 
=1.53*1/1.204*1.1018 
= 1.4 

From IFT corresponding to A e /A*=1.4, we have 
M e =0.47 

P 02 /P e =l/.859=1.163 
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P e =l/1.163*.917*l 
P e =.788 atm. 

But it is given that P e =.75 atma . , hence we to move the shock a bit downstream. 
STEP 4: 

Let us assume shock wave is located at A s /A t = 1.301 

We have M 1 =1.66,from NST corresponding to Mi=1.66, we have 

M 2 =.6512 and P 02 /P 0 i=.872 

From IFT corresponding to M 2 =.6512, we have A s /A =1.1356 

.A.g/.A. — .A.g/.A.^'^.A.^/.A.g'^.A.g/.A. 
=1.53*1/1.301*1.1356 
=1.335 

From IFT corresponding to A e /A*= 1.335, we have 
M e =0.5 

Po 2 /P e =l/.843=1.1862 
P e =l/1. 1862*872*1 
P e =.735atm 

Again the pressure ratio does not exactly match our given exit pressure. Now we can interpolate between the two assumed 
values. So 

A s /A t =1.301-(1.301-1.204)*(.75-.735)/(.788-.735) 
A s /A t = 1.274 
STEP 5: 

Let us assume that shock wave is located at A s /A t = 1.274 

From IFT corresponding to A s /A t = 1.274, we have Mi=1.63 

From NST corresponding to Mi=1.63,we have M 2 =.6595 and P 02 /P 0 i=.8838 

From IFT corresponding to M 2 =.6595,we have A s /A =1.1265 

.A.g/.A. — .A.g/.A.£ * .A.^/ .A. § * .A. s/.A. 
=1.53*1/1.274*1.1265 
=1.353 

From IFT corresponding to A e /A*=1.353, we have 
M e =0.49and P 02 /P e =l/.848=1.178 
P e =l/1.178*.8838*l 
=.750atm 

Thus the properties obtained for a C-D nozzle of area ratio A e /A t =1.53 at operating pressure ratio of 0.75 are 

Mi=1.63 

P 02 /Poi=.8838 

M e =0.49 

P 02 /P e =1.178 

P e =0.75 

A s /A t = 1.274 

III. RESULTS AND DISCUSSIONS 



TABLE. 1 FOR AREA RATIO 2.035 
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0.90 


0.84 


0.78 


0.69 


0.6 


0.5 




M e 


.312923 


.334810 


.359950 


.405523 


.464070 


.552165 




Po2/Pe 


1.070239 


1.080692 


1.093671 


1.119925 


1.159045 


1.230187 




P02/P0I 


.963215 


.907781 


.853063 


.772748 


.695427 


.615094 


hock 




1.380000 


1.566000 


1.707000 


1.889000 


2.055000 


2.229000 


(Py-Px)/Px 


1.055133 


1.694415 


2.232824 


2.996374 


3.760196 


4.629848 
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T 2 /T l 


1.241814 


1.364657 


1.463361 


1.599362 
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1.882576 


s acr 


P 2 /Pl 


2.055133 


2.694415 


3.232824 


3.996374 


4.760196 


5.629848 


ertie 


Pl/Pl 


1.654945 


1.974427 


2.209177 


2.498731 


2.747271 


2.990502 


Prop 


A s /A t 


1.104193 


1.223588 


1.344237 


1.541706 


1.767480 


2.057396 
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Fig. 3.1 variation of shock location with area ratio 
TABLE 2. Comparison between analytical and hit & trial method 



Properties 


Analytical method 


Hit and trial method 


Exit mach number(M e ) 


.3812 


.38 


Stagnation pressure ratio (P02/P01) 


.6632 


.6628 


Inlet mach numbe^Mi) 


2.1240 


2.12 


Temperature ratio(T 2 /Ti) 


1.7908 


1.787 


Pressure ratio(P 2 /Pi) 


5.0966 


5.077 


Shock location(A s /A t ) 


1.8755 


1.869 


Shock strength(P v -P x )/P x 


4.096 


4.077 



The above values have been calculated for properties across Normal shock waves, where working fluid consider 
here is air with specific heat ratio y=lA for a nozzle geometry of area ratio of 2.494 and the operating pressure ratio of 0.60 
which has been chosen arbitrarily in between the critical pressure ratios to compare the result obtained by both the methods 
i.e. "Analytical method" and "Hit& trial method". 

IV. CONCLUSIONS 

1) Shock strength increases significantly by increasing (Mi), whereas shock location do not vary much with (Mi). 

2) For constant operating pressure ratio, shock strength shows significant variation for smaller area ratio(A e /A t ). 

3) For a constant area ratio (A e /A t ), shock location varies significantly for higher pressure ratios as compared to lower 
pressure ratios. 

4) Shock location moves towards exit, by decreasing (P e /P 0 )- 
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